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Abstract: At a meson melting transition in holographic QCD, a weak turbulence 
of mesons was found with critical embeddings of probe D-branes in gravity duals. 
The turbulent mesons have a power-law energy distribution Sn oc (cUn)" where cOn is 
the mass of the n-th excited resonance of the meson tower. In this paper, we find 
that the turbulence power a is universal, irrespective of how the transition is driven, 
by numerically calculating the power in various static brane setups at criticality. 
We also hnd that the power a depends only on the cone dimensions of the probe 
D-branes. 
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1. Turbulence on branes in holography 

AdS instability conjecture by Bizon and Rostworowski [1] has attracted much atten¬ 
tion on intrinsic turbulent nature of AdS spacetime. The importance of addressing 
the stability question of generic AdS spacetimes is obvious from the viewpoint of the 
renowned AdS/CFT correspondence [2, 3, 4]. 

The issue of the AdS instability brings about various fruitful discussions on a 
possible universality of the instability (see for example [5]-[39]), while partially relies 
on details of numerical simulations. On the other hand, probe D-branes in AdS 
spacetimes, which are quite popular in AdS/CFT correspondence for introduction of 
quarks [40, 41, 42, 43, 44] and on which AdS spacetimes are generically induced, may 
provide an easier toy model for a similar turbulent instability. In fact, in Refs. [45, 46], 
turbulent spectra of meson energy are observed, which is quite reminiscent of the 
gravity instability in AdS spacetimes. The scalar and vector helds on the havor 
D7-brane in AdS^ x geometry correspond to mesons m. M = 2 supersymmetric 
QCD a la AdS/CFT, and a decomposition by the discrete momentum along the 
holographic (AdS radial) direction dehnes the energy distribution for the n-th level 
of meson excitation. It was found in [45] [46] that the energy distribution at a critical 
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brane embedding forming a conical shape with the help of background electric held 
exhibits 


En OC {Ujn) ^ ■ (1-1) 

Here Sn is the energy of the n-th resonance of the mesons, and ujn is the mass of 
the n-th resonant meson. Since the mass Un is nothing but the discrete momentum 
along the holographic direction, the critical behavior (1.1) is a weak turbulence. 

The interesting points in [45] [46] about the brane turbulence are that the turbu¬ 
lence (1.1) appears not only in time-dependent simulations of applied electric helds 
[47, 48] but also for a static brane conhguration under a static external electric held, 
which is called critical embedding [49, 50] (see [51] and [52] for the critical embedding 
at thermal phase transition). This is the reason why the brane turbulence may be 
tractable in a easier manner. 

In this paper, we analyze the weak turbulence on the probe D-branes in AdS 
geometries, and hnd a universality of the turbulent behavior, and also a universality 
of the exponent in the relation 


En OC {UnT ■ (1.2) 

Turning on electric helds and magnetic helds, and also temperature, we hnd that 
this power a does not change at criticality. It is dependent only on dimensions of 
the probe brane. 

Let us describe a little more on what is the criticality where the turbulent behav¬ 
ior (1.2) appears. The quark-hadron transition in strongly coupled gauge theories 
including QCD is one of the famous transitions, and can be realized either by turn¬ 
ing on a nonzero temperature, or by putting a strong electric held, or by some other 
methods. The phase boundary separates the hadron (“conhning”) phase and the 
quark (“deconhning”) phase. In large gauge theories where N^. is the number of 
colors, the phase transition can be hrst order. In the AdS/CFT analysis in a probe 
approximation of a havor brane, there appears a hrst order phase transition called 
meson melting [51]. When the temperature, or the external electric held applied to 
the quarks, is not large, the theory has a discrete spectrum of mesons. On the other 
hand, when it is large, the meson spectrum becomes continuous, and mesons melt. 
The transition shares properties with the quark-hadron phase transition in QCD 
(although generically gluons are always deconhned in this meson melting transition). 

The meson melting transition has been known to be accompanied by a conic 
D-brane conhguration (the critical embedding). It appears at the phase boundary, 
and should characterize the meson melting transition. The conic shape of the probe 
D-brane can be rephrased in terms of condensation of mesons, as the position of 
the D-brane in the bulk gravity geometry is probed by the expectation value of 
the mesons. In fact, for the case of a D7-brane in static electric helds, the critical 
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embedding shows [45, 46] the turbulent condensation of mesons characterized by 
(1.1). The relation (1.1) was also kept for dynamical simulation for time-dependent 
electric held backgrounds or time-dependent quark mass backgrounds [45, 46]. 

In this paper, we examine the relation (1.2) for various static backgrounds pro¬ 
viding critical embedding of probe D-branes. In Section 2, we consider generic electric 
and magnetic held on probe D7-branes in AdS, and also a temperature, and hnd a 
universal power 


a = —5 (D7 in AdS^). (1.3) 

In Section 3, we examine the case of a probe D5-brane in similar external helds, and 
hnd a universal power 


a = —4 (D5 in AdS^). (1.4) 

From these results and also some simple analysis of conic branes of 1 and 2 dimensions 
in hat space, in Section 4 we make a conjecture for generic brane turbulence 

Sn cx (a;„)-'^---^, (1.5) 

where dcone is the number of the spatial dimensions of the cone of the brane. 

2. D3-D7 system and universal turbulence 

In this section, we analyze the meson turbulence in the D3-D7 system with an elec¬ 
tromagnetic held or at a hnite temperature, by using the AdS/CFT correspondence. 
The turbulent behavior was found in this probe D7-brane system in [45], and here 
we explore diherent parameter space as external helds, to hnd that the turbulence 
power spectrum found in [45] is universal: a = —5. 

First, we review the D3-D7 system and the scalar meson mass spectrum. Next, 
we analyze a relation between the energy of the meson Sn and the meson mass ujn 
by changing the electromagnetic held and the temperature. We will hnd that the 
universal power-law of the energy of the higher states of the meson resonances goes 
as Sn oc when the probe D-brane is at the critical embedding. 

2.1 Review of the J\f = 2 supersymmetric QCD in AdS/CFT 

First, to hx our notation, we briehy review the ehective action of mesons for the 
four-dimensional M = 2 supersymmetric QCD given by a gravity dual of the D3-D7 
system [40]. In order to use the AdS/CFT correspondence, we consider the gauge 
group U{Nc) and one-havor quark sector with a large Nc limit Nc ^ oo and the 
strong coupling limit A = NcQym The brane conhguration of the D3-D7 sys¬ 

tem is shown in the following table: 
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In the gravity dual, the effective action of the meson tower is the flavor D7-brane 
action in the AdSs x background geometry, 

^ = fO ^6 I /8 / det{gab[w] + 27iP^Fab), (2.1) 

[ZTl) J 

ds^ = -^'qavdx^dx'' ^^ + p^dVt\ + div^ + , ( 2 . 2 ) 

where = p‘^ + w^+w‘^, Fab = da^ — db^a-, and R = {2XY^Hs is the AdSs radius. We 
treat the flavor D7-brane as a probe, so it does not affect the AdSs x background 
geometry, which is justified at Nc S> 1. For our purpose we can assume that the 
fields of the action depend only on x^ and p such that w = w{x^,p), w = w{x^,p) 
and Aa = Aa{x^,p). These fields correspond to the scalar mesons and the vector 
mesons, after they are decomposed as shown below. 

In order to solve the equation of the motion for (2.1), a boundary condition 

w{x^, p = oc) = R^m, (2.3) 

is necessary (while we set w{x^,p) = 0 by using a rotation symmetry of the {w,w)- 
plane). In the AdS/CFT correspondence, m is related to the quark mass nig as 
nig = (A/27r^)^/^m. The static solution of the equation of the motion for (2.1) with 
(2.3) is 


w{x^,p) = R^m, Aa{x^,p) = 0, (2.4) 

which is a straight D7-brane located at w{x^,p) = R^m. 

From now, we consider a scalar fluctuation of w{x^,p) around the solution (2.4) 
[41]. We define the scalar fluctuation field as x{t, p) = R~‘^w — m where for simplicity 
we assume that it does not depend on x^,x‘^,x^. Substituting it to (2.1), we obtain 
an action of y as 

f r°° P^ 

S= dtd?x / dp—r^ -——^ 

J io 2(p2 + i?4^2)2 

up to an overall factor. 

In order to determine the meson mass spectrum, 
motion for y to the first order, 

(9^ (p^ + d p^ d 

dF p^ dp dp 


we solve the equation of the 

y = 0. (2.6) 


{dtxf 


{p^ + R^rn^^^ 


■{dpxf 


+ 0{x^), (2.5) 
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(2.7) 


A normalizable solution of (2.6) is 

OO 

X = Re [Cn exp [iu:nt] Cn(p)], 

n=0 

en{p) = V2(2n + 3)(n + l)(n + 2) ( ^ 2 ^^ 4^2 
cJn = 2^J {n + 1 )(? 7 , + 2)m, 


n+l 

F{—n — 1, —n, 2; 


R^m'^ 




( 2 . 8 ) 

(2.9) 


where F is the Gaussian hypergeometric function. The eigen value ojn corresponds 
to the meson mass at the resonance level n. 

For the eigenmode expansion, we dehne the following inner product in the p- 
space, 

(/>«)-/ 'ip + ( 2 . 10 ) 

Under the inner product (2.10), e„(p) forms an orthonormal basis as 

(^n) Cm) ^nm- (2-11) 


By using this basis, we expand the scalar fluctuation held y as 

OO 

X = ^Cn{t)en{p). ( 2 . 12 ) 

n=0 

The expansion coefficient c„(t) corresponds to the meson held at the resonance level 
n. 

Substituting (2.12) to (2.5), we obtain the ehective action of the meson held 
Cn{t) as 


S = 




n=0 

where we omit total derivative terms, 
meson based on this action as 


— + interaction, (2-13) 

We dehne a linearlized energy of the n-th 


and the linearized total energy as 



(2.14) 


OO 


S — y ^ £n- 
n=0 


(2.15) 


Now, we are ready for the analysis of the meson turbulence. We shall see how the 
total energy is distributed to the n-th energies Irrespective of the external helds, 
the power law cx uj~^ at the criticality will be found below. 
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2.2 Turbulence with an electromagnetic field 

The turbulent meson condensation found in Ref. [46, 45] at static D7-brane config¬ 
uration is about a critical electric field, and the power a appearing in the scaling 
behavior Sn oc (a;„)" is found to be a = —5. Here in this section, we include a 
constant magnetic field in addition, and will find that again the power a is —5 at 
criticality. 

We introduce a generic constant electromagnetic field in the four-dimensional 
x^-spacetime and calculate the power-law scaling of the energy of the highly excited 
meson states at the critical embedding. We turn on generic constant electromagnetic 
fields as 


ri 2 — —-^21 = -^23 — 

and again consider a symmetric D 7-bra 
' ' magnetic field, we can calculate 

(p2+M;2)2^l I 


; the flavor 



\B\ 


0 



(p^ -f ra¬ 
the brane embedding w{p) is 


The Euler-Lagrange equation for 


III 

O 

CO 

(2.16) 

-Bi 3 = By, 

(2,17) 

- w{p), w = 0. 

With this 

action as 


R8(2p/2). 

(p2 -i- w'^y 

liE-Sy. 


(2.18) 


(p2+^2)2 


.,/2 




,,/2 


2R^(27r;g)^w / I pi: 

p2+^2)3 IRI 


I I 


y (p2+u!2)2 


a/TT” w'‘^ 


\B\' 


(2.19) 


We calculate the energy distribution of the meson resonances by the following 
steps: 


1. Solve (2.19) with (2.3) and calculate x{p) numerically. 

2. By using (2.11) and (2.12), calculate Cn from x(p). 

3. By using the definition (2.14), calculate Sn from c„. 

For numerical calculations, we set i? = m = 1, and we normalize Si = 27r/^i7j and 
Bi = 2'Kl‘j.Bi. In particular, we are interested in a relation between and Un at the 
critical embedding of the probe D7-brane. 

In Fig. 1, Fig. 2 and Fig. 3, we plot numerical solutions of w{p), for given Lorentz- 
invariant combinations |Tp — \B\^ and £ ■ B. We find the critical embedding (conical 
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D7-brane configurations) which are of our interest. In each hgure, the red curves are 
the critical embeddings. 

Fig. 1, Fig. 2 and Fig. 3 correspond to different values of \£\^ — respectively. 
In each hgure, for those hxed values of \£\‘^ — we vary £ ■ B and hnd the critical 
embedding (red curves in each hgure). We dehne {£ ■ B)ct as the value oi £ ■ B which 
gives the conical solutions.^ 

w{p) w{p) 



Figure 1: The shape of the D7-brane with \£\‘^ — \B\'^ = 0. The curves correspond 
respectively to £ ■ Bl{£ ■ B)cr = 0.98, 0.99, 1, 1.001, 1.001, 1 from top to bottom and 
{£ ■ B)ct ~ 0.45793. Right figure is an enlarged view around p = 0. 


w{p) w{p) 



Figure 2: The shape of the D7-brane with |Tp — = —0.25. The curves correspond 

respectively to £ ■ Bl{£ ■ B)cr = 0.98, 0.99, 1, 1.001, 1.001, 1 from top to bottom and 
{£ ■ B)ct ~ 0.63647. Right figure is an enlarged view around p = 0. 


^Note that there are solutions whose absolute value oi £ ■ B is larger than that of {£ ■ B)^. 
Furthermore, for given electromagnetic fields, there may be several brane solutions. It is expected 
that those solutions are continuously connected to the critical embedding in the parameter space, 
with a fractal-like structure (see some discussions in Ref. [49]). Here we are just interested in the 
critical embedding. 
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w{p) w{p) 



Figure 3: The shape of the D7-brane with \E\^ — = +0.25. The curves correspond 

respectively to £ ■ i3/(T • B)cr = 0.98, 0.99, 1, 1.003, 1.003, 1 from top to bottom and 
{£ ■ 13)ct ~ 0.21273. Right figure is an enlarged view around p = 0. 


Finally we present the power-law of the energy distribntion. Fig. 4, Fig. 5 and 
Fig. 6 are plots of log(e:„/£o) as a function of log {un/m). We notice that the red 
dotted lines (which are for the critical embeddings) are described well by a linear ht, 
in particular at larger regions, in each hgure, while for other smooth embeddings 
the ht cannot be linear and the energy decrease rapidly for larger The linear ht 
for the red dots means the power-law of the energy distribution which we expected. 
Black lines are htting lines of the red dots at large a;„. From the slope of the black 
line, we obtain the power-law scaling of the energy distribution as 

oc (2.20) 

Since the power appearing in the ht does not depend on |Tp — of our choices, 
the power is universal, and coincides with a = —5 given in Ref. [46, 45]. 


2.3 Turbulence at a finite temperatnre 

In the previous subsection, we have found the universal power law (2.20). Here, 
instead, let us turn on a temperature rather than the electromagnetic held to hnd 
the universality of the power law. We introduce a nonzero temperature following 
[54], and calculate the power-law scaling of the energy distribution at the critical 
embedding. 

A nonzero temperature can be easily introduced by replacing the background 
AdS geometry by the AdS black hole metric as 


= ^[-f{r)dt‘^ + dx"^] + ^ 


dr"^ 

J{r] 




( 2 . 21 ) 


peculiar feature of Fig. 5 is that ei is larger than erg. This is related with the shape of 
the D7-brane in Fig. 2 where dpw{p) is negative at large p, which shows the repulsion due to the 
world volume magnetic field. 












Figure 4: Plot of log(e„/eo) as a function of \og{ojn/m) with \£\’^ — \B\'^ = 0. The color 
of dotted lines corresponds with that of Fig. I’s curves. Black line is a fitting line of red 
dots at large ujn- 



Figure 5: Plot of log(e„/eo) as a function of log(w„/m) with |Tp — |i3p = —0.25. The 
color of dotted lines corresponds with that of Fig. 2’s curves. Black line is a fitting line of 
red dots at large Un- 


where /(r) = 1 — {th/tY and the location of the horizon rn is related to the temper¬ 
ature T by the AdS/CFT dictionary as T = We introduce a new coordinate 

u as 


2u = r + — rjj, 

with which the AdS black hole metric (2.21) is written by 


2 r 


. 2 “ 
* =ip 


u 


Ru) 

RyR 

Rn) 


dR -|- f{u)dx^ 


dtR + f{u)dx‘^ 


+ + u^d£ll] 




where 




( 2 . 22 ) 


-I- -[dv^ + v^dVtl +dw^ + dw^]^ (2.23) 


(2.24) 
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Figure 6: Plot of log(e„/£o) as a function of log(t(j„/m) with \£\^ — \B\‘^ = +0.25. The 
color of dotted lines corresponds with that of Fig. 3’s curves. Black line is a fitting line of 
red dots at large Un- 


As in Sec. 2.2, we use the ansatz 


w = w{v), tc = Aq = 0, 


and a boundary condition 


w{v = oc) = R^m. 


(2.25) 


(2.26) 


Substituting our ansatz (2.23) to (2.1), we obtain the flavor D7-brane action as 


duu'* 1 - 


S' = y dt(rx 

The Euler-Lagrange equation for w{v) is 

d„w 


H 


16m® 


V^l + (dyw)'^. 


(2.27) 


0 = a, 


' H 


16(^2+it;2)4y + {dyWy 




V wr 


H 


2(^2 + 


Vl + {d,wY. (2.28) 


Following the strategy of [45, 46] and that in Sec. 2.2, we calculate w{v) and from 
(2.28) with (2.26) numerically. For the numerical simulations we set R = m = 1 as 
before. 

Fig. 7 describes the shape of the D7-brane with the finite temperature effect. To 
find a critical embedding which is a conical shape of the flavor D7-branes, we vary 
the value of th while fixing the quark mass m = 1. The red conical solution in Fig 7 
is the one which we want. We define (rj^)cr as the value of th which provides the 
conical solution. In the figure we also plot other solutions of Minkowski embedding 
for comparison.® 

^Note that lowering the tip of the D7-brane does not directly corresponds to raising the tem¬ 
perature. For a certain region of values of there exists several solutions for a given value of ru- 
Furthermore, there are solutions whose is larger than (r//)cr. The situation is similar to that in 
the previous section. The existence of such solutions was discussed in Ref. [54]. 
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w{v) w{v) 



Figure 7: The shape of the D7-brane with the finite temperature effect. The curves 
correspond respectively to rH/{TH)cr = 0.98, 0.99, 1, 1.0005, 1.0005, 1 from top to bottom 
and (r//)cr ~ 1.0907. Right figure is an enlarged view around u = 0. 


Fig. 8 is a plot of log (e„/£o) as a function of log (un/m), for our probe D7-brane 
at finite temperature. The notation of the colors follow that of the previous section; 
the red dotted line (corresponding to the critical embedding) is found to obey a 
linear £t at large Un region. It means a power-law of the energy distribution for 
highly excited meson states at the critical embedding. The black line is a linear fit 
of the red dotted points at large Un, and we hnd a power-law scaling again, 

oc (2.29) 



Figure 8: Plot of log (en/so) as a function of log (ta^/m) with the finite temperature effect. 
The colors of dotted lines correspond to those of Fig. 7’s curves, respectively. The black 
line is a linear fit of the red dots at large Un- 


Summarizing the above, we have found the power-laws of the energy distribution 
for highly excited scalar meson states as (2.20) and (2.29). The powers are very close 
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to a = —5. Thus we find the universality,"^ 

En oc . (2.30) 

The power scaling is universal at critical embedding, and does not depend on how 
the critical embedding is realized: either by various electromagnetic fields or a finite 
temperature. 

3. D3-D5 system and universal turbulence 

In this section, we analyze the turbulent behavior of the mesons at excited states in 
D3-D5 brane system by using the AdS/CFT correspondence. The purpose is to find 
(1) whether we find the power-law behavior, and (2) whether the power is universal. 

We study two ways to obtain the critical embedding of the probe D5-brane, as 
in the previous section. After reviewing the energy of the meson at zero temperature 
and with no electromagnetic field, we introduce a constant electric field, and study 
how the relationship between the energy spectrum and the mass of the mesons 
is by changing the background values and accordingly the shape of the D5-brane. In 
addition, we are interested also in the turbulent power law for the meson states at a 
finite temperature, and we examine the universality. 

In this section we will find the energy distribution as En oc at critical 

embeddings of the D5-brane, either for a constant electric field or a finite temperature 
system. The power is shared by both the situations, so is for the D3-D5 system. 
However, the power itself is —4 which is different from the value of the power, —5, 
for the D3-D7 system. 

3.1 Review of the M = 2 supersymmertic defect gauge theory in AdS/CFT 

In this subsection, we review the derivation [53] of the spectrum of “mesons” from 
the fluctuation of a probe D5-brane at zero temperature and with no background 
electric field, by using the AdS/CFT correspondence. We are interested in the fol¬ 
lowing D3-D5 brane configuration: 
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The brane configuration preserves M = 2 supersymmetries in total. The flavor probe 
D5-brane extends along the directions x®, x®, so it shares only x®, x^, x^ 

directions with the gauge Nc. D3-branes. Thus, the gluon A/” = 4 multiplets live in 
the (3-1-1) dimensional spacetime while the quark hypermultiplets (and resultantly 

^We expect that the numerical power-law scaling factor will approach —5 if we extend our 
analysis to much higher states. 
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the mesons as their bound states) live only in a (2+1) dimensional domain wall, 
which is a defect. 

It is known that meson states can be analyzed holographically by the D3-D5 
brane system. The scalar meson held corresponds to a huctuation of the probe D5- 
brane along the directions due to the AdS/CFT correspondence. For the 

huctuation, the Laplace equation is classically solved by a Gauss hypergeometoric 
function [53], as we will see below. 

The AdSs x background metric is 

ds‘^ = -^'q^^dx^dx" H—y [dp^ + p^dVt^ + dul + du'l + dwgj. (3.1) 

The 5-sphere together with the AdS radial direction r is decided into a combi¬ 
nation of p and 2-sphere G 2 and the remaining three directions x'^,x^,x^ which 
are parametrized by a; 4 ,a; 5 ,a; 6 - So the relation among these parameters is = 
+a;| + a;| Tcug. Since (a; 4 , cug, cjg) has a rotation symmetry, we may hx cug = cug = 0 
for simplicity. 

The D5-brane action is a Dirac-Born-Infeld(DBI) action given by 

^DBi = -rs J e-V-det(P[p]ab + 27r/2F,,), (3.2) 

where rg is the D5-brane tension and is dehned by rg = l/(27r)^Ps/®. Qs is the string 
coupling, (j) is the dilation held which is set to zero as the background. 

The scalar mesons are eigen modes which are liberalized solutions for huctuations 
of u}i{x^, p) of the D5-brane DBI action to which the background metric is substituted. 
Denoting the index i running i = 0,1, 2 as the D5-brane does not extend along the 
direction x^, we impose a boundary condition at the AdS boundary as 

a ;4 (P, p = cxo) = i?^m, (3.3) 

where m is related to the quark mass rrig as nig = (A/27r^)^Am due to the AdS/CFT 
dictionary. A static classical solution of the DBI equation of motion is 

U4{x\p) = R^m. (3.4) 


The solution roughly measures the distance between the D3-branes and the D5-brane. 

Next, we consider the huctuation y around the static solution i?^m, dehned by 
X = — m. We assume for simplicity that y is independent of the coordinates 

x^ and x"^. The action for y obtained by just expanding (3.2) to the quadratic order 

is 


S 



p^R^m 

2(p2 + R^m?Y 



(p2 + R^m'^f 



+ (^(x"), 


(3.5) 
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where the irrelevant overall factor is neglected. From (3.5), we derive the equation 
of motion as 


d'^ {p^ + R^7ii?Y d p^m d 
dt^ p'^R^m dp dp 


X = 0, 


Its normalizable solution is 


X = ^ Re [Cn exp[innt]Enip )], 


n=0 


The basis functions are given by 


(3.6) 


(3.7) 


4(n+l)/ R‘m^ \ , , 

where F is the Gaussian hypergeometric function. The mass of the level n resonance 
meson is 


= 2^(l/2 + n) (3/2 + n) m. (3.9) 

For our later purpose we dehne the inner product in the p-space as 

roo 2 p2™ 

‘''’(yTR%(T^’(^')G(p). (3.10) 

under which we have the orthonormality condition 

{En,Em)=5nm- (S-H) 

Using the basis (3.8), we expand the fluctuation of the scalar held as 

OO 

X = ^Cn{t)En{p), (3.12) 

n=Q 

then with this we can dehne the linearlized meson energy at level n, 

^n = ^icl + nlcl). (3.13) 

The linearlized total energy is given by the formula (2.15). All of these are analogous 
to what we have used in the case of the D3-D7 brane system in Sec. 2. 

3.2 Turbulence with an electric field 

We are ready for studying the turbulent meson behavior. The turbulence should 
show up in the meson energy distribution which is dehned as (3.13). Note that the 
energy spectrum (3.13) is dehned with no background electric held. Once the electric 


14 










field is turned on, the meson spectrum changes compared to that with no electric 
held, since the probe D5-brane is affected by the electric held. It is well-known 
that the shape of the probe D-brane is deformed by the electric held [49]. We are 
interested in how the shape of the D5-brane changes and how the meson spectrum 
is ahected accordingly. 

To look at how the D5-brane shape is deformed, we solve the equation of motion 
from the DBI action (3.2). AdSs x background is given by (3.1). From a simple 
consideration of rotational symmetry, we can put cus, = cue = 0 while depends on 
the radial direction p as = L{p). Then, the induced metric on the D5-brane is 
given by 

+ ^ijdx'dx^) + ^[{1 + {dpVf}dp^ + ( 3 - 14 ) 

where i,j = 1, 2. Then we turn on a constant electromagnetic held. In (2-1-1) dimen¬ 
sions, Lorentz transformation can bring us to a frame on which only Foi component 
is nonzero. Denoting E = Fqi, we obtain the action with the electric held as® 


5'nRT = -dTTT^ / (fx 


dp pyi+1 - (3.15) 


The equation of motion is obtained with a redehnition S = 2'iil‘j.E as 

(p5WF 2R^8^p^L^l + {dpLf 


dn 


p^dpLJl - 


V^l + {d)pLy 




= 0 . 


(3.16) 


We solve this equation numerically to obtain the shape of the probe D5-brane for a 
given value of £. The results of the numerical calculations, with R = 1 and m = 1, 
are shown in Fig. 9. the shape of the D5-brane changes and the D5-brane has a cusp 
as the static electric held increases. With various chosen background electric held 
the D5-brane changes its shape, and we hnd a conical D5-brane. It is the critical 
embedding, for which we expect the turbulent behavior of mesons. 

Basically, as in the case of the D7-brane probe, for small £ the D5-brane is at the 
Minkowski embedding (which has a smooth shape), while for a large enough £ the 
D5-brane is at a black hole embedding where the induced metric on the D5-brane 
has an ehective horizon. The critical embedding is in between those embeddings, 
and the D5-brane has a cone at the tip. The value of £ with which the D5-brane 
can be in the critical embedding is called £’„, and in Fig. 9 we plot various curves 
for various £/£cv^ 

®The Chern-Simons terms in the D-brane effective actions does not contribute in our analysis. 
For D4-D8 setup, they may contribute. See a discussion in [58]. 

®Note again that the location of the tip of the D5-brane is not a monotonic function of EjEcp- 
In some region of the value of f/fcrj the D5-brane profile is found to be not unique. A fractal-like 
structure emerges, as in the case of the D3-D7 system. 
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Figure 9: The shape of the D5-brane in the AdS background. Each curve corresponds, 
from top to bottom, to S/Sct = 0.9,0.99,1,1.017,1 respectively. The D5-brane can have a 
cusp at p = 0 for T = Tcr(— 0.437) (the red curve). 



Figure 10: The power law of the re-th meson mass including a static electric field. The 
vertical axis is the logarithm plots of the n-th meson energy divided by the lowest meson 
energy. The transverse axis describes the logarithm plots of the n-th meson mass. 


From these numerical shape of the D5-branes in Fig. 9, we can calculate the 
meson energy spectrum, through the definitions given in the previous subsection. As 
the shape of the D5-brane changes by varying S/Sa, the decomposed meson energy 
spectrum changes. Our result for the energy distribution Sn is presented in Fig. 10. 

The energy distribution Fig. 10 shows that the critical embedding (the red dots) 
is distinctively different from other Minkowski embeddings. The red dots can be 
linearly £t as 

cx (On)-'-"*" (3.17) 

which is a power law. This is nothing but a weak turbulence, and is very similar 
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to what we has been known in [45, 46] and what we found in the previous section. 
For the other Minkowski embeddings, there is a signihcant reduction of energy for 
higher excited mesons (large n region). So, the power law appears only at the critical 
embedding. The cusp of the D5-brane seems to be responsible for the power law of the 
meson energy distribution. We can conclude that the weak turbulence of the excited 
mesons is caused by the cusp of the D5-brane, which is realized and accompanied 
with the phase transition. 


3.3 Turbulence at a finite temperatnre 

In this subsection, we consider the energy spectrum of the highly excited meson 
states in a hnite temperature D3-D5 system in the absence of the electric field, to 
study the universality of the power law (3.17). The basis of the analysis was given 
in [54] as we have reviewed. In the previous subsection, we studied the turbulent 
meson spectrum in the electric field. When the probe D5-brane has a cusp, the meson 
energy distribution obeys a power law: Sn oc It is of interest if the power 

law is universal or not. So, here we introduce a temperature to the AdS background, 
and examine the power law. 

In AdS/CFT correspondence, the temperature is introduced by replacing the 
background AdS geometry (3.1) by an AdS black hole metric. 


ds^ = ^[-f{r)dt^ + dx^] + ^ 


dr'^ 

J{r] 


+ r^dVLi 


(3.18) 


where dx'^ = dx\ + dx 2 + dx‘l, and the function of /(r) is dehned by /(r) = 1 —(rH/r)"^. 
The location of the black hole horizon rn is related to the temperature as T = r-a/nR^ 
due to the AdS/CFT dictionary. Using a new coordinate 2v? = + \/r^ — r^, the 

AdS black hole metric is written as 


ds'^ 



/(^)^ 7.2 I fr \j^2 

—dt +f{u)dx 

L f{u) 



[dn T V dVl2 T dco^ T dcug T dcug] , 


(3.19) 


where 

+ (3.20) 

The radial coordinates u and v parameterize the radii of and respectively, and 
they are related as + 0)4 + 0)5 +a)g. As before, using the rotational symmetry, 

we can restrict ourselves to 0)5 = a)g = 0 and 0)4 = L{v). Then, the induced metric 
on the D5-brane is given by 


7 2 “ 

* =5? 


2 r 


fiu) 


dt^ + f{u)dx'^ 




{1 + {dyL)'^}dv'^ + v'^dQl 


(3.21) 
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Figure 11: The shape of the D5-brane in the AdS black hole background. Each curve is 
with rn/rcr = 0.9,0.99,0.99,1.01,1 respectively from top to bottom. The D5-brane of the 
red curve has a cusp at p = 0, which we call a critical embedding, for th = rcr(— 0.443). 
The other curves are Minkowski embeddings. 


Next, we determine the shape L(v) of the probe D5-brane in this finite temper¬ 
ature system. The 1-flavor DBI action is defined by (3.2). By using the induced 
metric (3.21), the DBI action becomes 


5'dBI — —dvTTs 




4(u2 + L2)2 


1 + (a„L)2. 

(3.22) 


The Euler-Lagrange equation of the D5-brane is obtained as 


v'^dyL (4(u^ -|- I/^)^ — ’"h ) V 4('y^ + + f' 


(3.23) 


(u2 + L2)3Dl + (a„L) 


24 ^ 2 ^ 4(u2 + l 2)2 + 34 ./i + (a,L)2 


= 0 . 


(u2 + L2)4.i/4 (u2 _|_ + r 


H 


We numerically calculate classical solutions L{v) of the equation of motion. We chose 
the convention R = 1 and m = 1 which is the same as before. The theory has only 
two scales, the quark mass and the temperature. So the theory is determined only by 
the ratio of those. For fixed m = 1, we vary th and obtain various D5-brane shapes, 
as shown in Fig. 11. 

As the temperature changes, the amount of the gravity which the D5-brane feels 
changes, since the location of the horizon approaches the D5-brane. The red curve 
in Fig. 11 has a cusp of the D5-brane, and shows that the D5-brane is at the critical 
embedding. 
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We expect that a turbulent behavior of meson excited states causes the cusp of 
the D5-brane at rn/rcr = 1 shown by the red curve in Fig. 11. So let us present the 
results of the meson energy spectrum for the higher excitation modes: see Fig. 12. 
As expected, the red dots which are the energy distribution of meson excitations 
corresponding to the cusp D5-brane (the red curve in Fig. 11) has a linear behavior. 
A linear £t of the red dots shows 


oc (3.24) 

Again, we have found that the critical embedding of the D5-brane shows a turbulent 
behavior of mesons. For other curves with the Minkowski embedding, we have no 
power law: the energy deposit at the higher meson resonance decreases rapidly for 
large n. So the weak turbulence, the power law behavior, is unique to the critical 
embedding of the D5-brane. 



Figure 12: The power law of the energy distribution for the n-th meson resonance, 
in the finite temperature system. The vertical axis is the logarithm of the n-th meson 
energy divided by the lowest meson energy eq. The horizontal axis is for the logarithm 
of the n-th meson mass fin- The colors of the dots correspond to that of the curves in the 
previous figure. 


In conclusion, the turbulent behavior of the higher excited modes of mesons is 
observed when the shape of D5-brane has a cusp in the static electric held or at 
the hnite temperature. The meson energy distribution obeys the power law as a 
function of the meson mass, a ~ —3.97 or /3 ~ —3.95 in the static electric held or 
at the hnite temperature system, respectively. Since these powers are quite close to 
—4, we conclude that the turbulent behavior is a universal phenomenon irrespective 
of how the transition is made. We expect that the exact power may be Sn oc (n^)""^, 
which may be conhrmed with numerical simulations with higher accuracy. 
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4. Universal turbulence and a conjecture 

In the previous sections, we found that the power associated with the level distribu¬ 
tion En of the energy density for the meson level n has an integer power: 

£„oc(mn)“, (4.1) 

where a = —5 for the D3/D7 system (treated in section 2) and a = —4 for the 
D3/D5 system (treated in section 3). Furthermore, the power in each case is found 
to be universal. The power does not change even though the external held to produce 
the phase transition is modihed — the temperature and the magnetic held. 

Therefore, it seems that the power a depends only on the dimensionality of the 
brane cone, and it does not depend on how the phase transition is driven. It is natural 
to make a generic conjecture on the power a\ The energy distribution for level n at 
the conical brane configuration for the phase transition is given by the power law as 
in (f-l), and the power is determined as 

a = —{dcone + 1) (4-2) 

where dcone is the dimension of the cone. For the case of the D3/D7 (D3/D5) 
system, dcone = 4 (dcone = 3). 

In this paper, we presented evidence for this conjecture for various situations 
associated with the D3/D7 and the D3/D5 holographic models, with various external 
helds, and all are consistent with the above conjecture. 

To strengthen the plausibility of the conjecture, in the following part of this 
section we provide simple examples of decomposition of a conic brane conhguration 
in flat space by eigen modes of harmonic functions. The examples are 1-dimensional 
and 2-dimensional cones, for simplicity. Both nontrivially can be worked out and are 
consistent with the conjectured relation (4.2). 

Ex. 1-dimensional cone 

Let us consider a one-dimensional cone in a flat space and calculate the power 
a when it is expanded in eigen modes. The one-dimensional cone is simply 
written as 

L{p) = 1 - -\p\ (4.3) 

71 

which is dehned on an interval —ti/ 2 < p < 7r/2, at whose boundary the height 
function L[p) is put to zero. We expand this cone by harmonic functions in 
one-dimensional flat space, i.e. Fourier modes. The eigen basis satisfying the 
Dirichlet boundary condition at the boundaries of the interval is 
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where n = 0,1, 2, • • ■ specifies the level of the modes. The magnitude of the 
eigenvalue of the Laplacian (9^ is just m\ = (2n + 1)^. The expansion is easily 
done as 

~ 2^/2 1 
L{p) = ^Cnen{p) where (4.5) 

So the meson energy for each level n is 

1 2 2 2^ 1 

= ^m^Cn -^^2n + 1)2 

So, for large n, we obtain the power 

En ^ n°‘, a = —2 . (4.7) 

From this simple example, we find that the possible power in the meson melting 
transition, if occurs in holography, would have a = —2, if the cone of the probe 
brane is just one-dimensional. 

Ex. 2-dimensional cone 

Let us proceed to an example in 2 dimensions, to support our conjecture. The 
example is again a cone in a fiat space. Let us consider a 2-dimensional cone 


Hp) = 1 - P 


(4.8) 


where p = a /+ (x^)^, and we consider 0 < p < 1. The boundary p = 1 is 
a circle, and we impose a Dirichlet boundary condition for harmonic functions 
on the disk. The harmonic functions in 2 dimensions are Bessel functions Jui'f'), 


dP Id 
dr‘i r dr 



(4.9) 


The subscript p corresponds to the angular mode, so our cone configuration 
should be expanded only with u = 0. The Dirichlet boundary condition says 
e„(p = 1) = 0, therefore, in terms of the Bessel functions, we have 

en(p) = I Mrnp) (4.10) 

where r„(n = 0,1, 2, 3, • ■ ■) label zeros of the Bessel function Jo(^)- The nor¬ 
malization is already fixed by using the following formula 

^ P {Joprnp)f dp=]^ (To(rn))^ • (4.11) 
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The eigenvalue of the eigen mode function en{p) is r^, thus the “meson mass” 
of level n is = Vn- 

The expansion of the cone (4.8) is given by 


oo „i 

L{p) = '^Cnen{p), Cn= p{l - p)en{p) ■ (4.12) 

n=0 

We are interested in only the large n behavior, so we use the following expres¬ 
sion for the large n asymptotic expansion of the Bessel functions, 


Jo(r) = Y— ^cos(r — 7r/4) -|- — sin(r — 7r/4) -|- 0{l/r^) 


The zeros of this function at large n are given by 


4 ? 7 , -|- 3 


-TT ■ 


(n > 1) 


4 ' (8n -|- 6)7r 

So, at the leading order in the large n expansion, we have 


[2 f {An + 3)7r tt 

en(p) ~ - - -p - - 


(4.13) 


(4.14) 


(4.15) 


Substituting this into (4.12), we obtain the large n behavior of the coefficient 

Cn. 1 




3^2 

47J-2 


5/2 


+ 0{n-^) . 


(4.16) 


The meson condensation c„ is found to scale as ~ n Therefore the energy 
distribution is obtained as 




^ 2 2 
— Til, C 
2 


9 -3 

16752 ^^ ’ 


(4.17) 


which shows the power 


CK = —3 

for the case of the 2-dimensional cone in a flat space. 


(4.18) 


These simple examples are consistent with our conjecture (4.2), so it is expected 
that the brane turbulence (4.1) with the power (4.2) is universal, not only in AdS-like 
geometries but also in a box of flat geometries. 
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5. Conclusions and discussions 


In this paper, we examined the universality of the power a which appears in the 
meson turbulence relation (1.2). The turbulent meson behavior was discovered in 
[45, 46] for a holographic system made by D3-D7 conhgurations which is dual to 
M = 2 supersymmetric QCD at large N^. and strong coupling limit. The power 
a = —5 was found in [45, 46] when the meson melting transition was driven by 
an external electric field, either static or dynamical. We found that the power a is 
universal, by calculating the power for the meson melting transition with magnetic 
fields in addition to the electric field, or with a finite temperature. For any transition 
we considered, we found a = —5 for the D3-D7 holographic system (Section 2), and 
a = —4 for the D3-D5 system (Section 3). 

Based on these findings with numerical calculations of a, we made a conjecture 
a = — (dcone + 1) where dcone is the number of dimensions of the cone formed by the 
probe D-brane. For the probe D7-brane (D5-brane), the cone dimension is dcone = 
4(3), respectively. 

Although we worked only with near horizon geometries produced by W D3- 
branes, we believe that other geometries with various other types of probe D-branes 
in holography may exhibit the same turbulent relation (1.2) with our universal power 
a = — (dcone + 1) determined solely by the dimensions dcone of the cone at the critical 
embedding. For example, it should be interesting to work out the power numerically 
for other geometries such as the ones made by Dp-branes {p ^ 3) or M-branes.^ 

The weak turbulence in holographic spacetime may be a quite generic phe¬ 
nomenon. In this paper we analyzed the turbulent energy distribution on probe 
D-branes, and found an intriguing universality. If the universality is applied to any 
probe sector of any holographic geometry, then we can make a generic statement 
that at “meson” melting transition in holographic field theories the meson energy 
distribution obeys a turbulent energy distribution (1.2) with a power universally de¬ 
termined by the dimensions of the probe conic brane. The dimensions apparently 
depends on the spatial dimensions of the internal space in the gravity dual. A radi¬ 
cal conjecture which follows from this argument is that the large QCD at strong 
coupling limit at a quark-hadron transition may have the meson turbulence with the 
power a = — 1, since QCD does not have any adjoint scalar field probing the internal 
space dimensions and thus it is expected to correspond to dcone = 0. 

The power a in the Kolmogorov-like power scaling of the turbulent relation (1.2) 
may be related to how easily the system is driven to the meson melting transition. 

^While writing this manuscript, we noticed that a recent paper [55] treats kinks on a fundamental 
string in a holographic setup. It may correspond to our dcone = 1 treated in Section 4, and the 
authors of [55] found numerically a ~ — 2 for a transverse circular quench, which is consistent with 
our formula. However, other quenches show different powers, so a further investigation for the 
universality may be necessary. 
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Once some external field which drives the system to the melted meson phase is 
turned on, the system becomes unstable and starts to evolve dynamically. Since the 
instability under strong electromagnetic field, as well as the time evolution, have 
been analyzed in holographic set-ups with D-brane effective actions (see for example 
[56, 57, 47, 58]), it would be interesting to investigate possible relations to those. 
Relatedly, it is important to ask what kind of physical properties in phase diagrams 
of strongly coupled gauge theories is characterized by the Kolmogorov-like power a 
of the meson turbulence. 
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